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Abstract 

The extended wings of the Call H&K lines provide excellent diagnostics 
of the temperature stratification of the photosphere of the Sun and of other 
cool stars, thanks to their LTE opacities and source functions and their large 
span in formation height. The aim of this study is to calibrate the usage of 
the H & K wings in one-dimensional interpretation of spatially averaged spec- 
tra and in deriving per-pixel stratifications from resolved spectra. I use multi- 
dimensional simulations of solar convection to synthesize the H & K wings, derive 
one-dimensional models from these wings as if they were observed, and compare 
r-j ■ the resulting models to the actual simulation input. I find that spatially-averaged 

O J models constructed from the synthesized wings generally match the simulation 

averages well, except for the deepest layers of the photosphere where large 
thermal inhomogencitics and Planck-function nonlinearity gives large errors. 
The larger the inhomogeneity, the larger the latter. The presence of strong net- 
work fields increases such inhomogeneity. For quiet photospheric conditions the 
temperature excesses reach about 200 K. One-dimensional stratification fits of 
discrete structures such as granulation and small-scale magnetic concentrations 
give satisfactory results with errors that are primarily due to steep temperature 
gradients and abrupt changes of temperature with depth. I conclude that stratifi- 
cation modeling using the H & K wings is a useful technique for the interpretation 
of solar high-resolution observations. 

Keywords: Models, photosphere - Granulation - Magnetic fields, photosphere 
- Spectral line, intensity and diagnostics 



1. Introduction 

The wings of the Ca n H & K lines at 3933 A and 3968 A span the solar photo- 
sphere in their formation which is relatively simple to interpret. They contain 
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suitable blend-free windows to sample intensities as well as blends to sample 
Dopplcrshifts at different heights. The span of their height sampling is extraor- 
dinary wide. |Ayres| (| 1989[) demonstrated that the opacity minimum at 4000 A 
and the large short- wavelength Planck sensitivity to temperature make the H & K 
background continuum respond as deep as the H~ opacity minimum near 1.6 /im. 
The H & K cores are chromospheric at line center and sample the shock-ridden 
int ernetwork "clapotisphere" as Hw and K?v grains near line center (e.g.. Rut- 
ten. I1995p . The extended wings sample all heights between these extremes, with 
clean single-peaked contribution functions because virtually all calcium resides 
in the Ca n ground state without sensitivity to Boltzmann and Saha population 
partitioning (Figure 6 of ILeenaarts et al\ I2006b|) and without sensitivity to 
NLTE population processes such as photon pumping or photon suction. The 
source functions of H & K also obey LTE-CRD until about 1 A from line cen- 
ter {e.g., |Shine, Milkey, and MiEdas) IT3751 |Ayres[ H3751 lUitenbroekl 11553)1 . 



IQwocki and A uer ( 19SUJ) emphasized the closeness of its source function to LTE 
which results from the relatively high probability of wing-photon destruction by 
incoherent scattering into the thermal line core. 

H & K wing modeling of solar and stellar photospheres was initiated during 
the 1970s at Boulder by J.L Linsky and collaborators, in particular T.R. Ayres 
QAyres, Linsky, and Shine) IT3741 [Ayres and Lmsky) fT975l |Ayres[ [T3751 11377)1 . 



Later, Rouppe van der Voort| ()2002p used the technique on a sunspot penumbra. 



The original application to solar faculae by Shine and Linsky ( 1974|) was followed 



by |Sheminova, Rutten, and Rouppe van der~ V oort (2005, henceforth Paper I) 



to derive models of individual magnetic concentrations ( "fluxtubes" ) from high- 
resolution H&K spectra. Can H wing modeling using an inversion SIR-code 
based on response functions was initiated by IBeck et al.\ 1)2008]) . It serves to 
interpret large sets of spatially resolved line profiles in an automated maimer. 
Recently, |Henriques and Kiselman] (|2009|) briefly described an application to 
high-resolution Call H filtergrams, including calibration using snapshots from 
numerical multi-dimensional simulations. 

Even more recent are the new results of lUitenbroek and Criscuohl (|201ip who 
demonstrate that the analysis of average spectra in terms of one-dimensional 
models is inherently unreliable because the inferred temperatures differ from 
the actual averaged temperatures. The average continuum of an atmosphere 
with temperature fluctuations is higher than the continuum of a one-dimensional 
atmosphere with the same average temperature stratification, due to the non- 
linearity of the Planck function with temperature. The interpretation of an 
average spectrum from a horizontally inhomogencous atmosphere in terms of a 
one-dimensional model therefore leads to an overestimate of the average temper- 
atures in the atmosphere. Since most of the information about solar and stellar 
atmospheres is obtained from unresolved spectra, it is important to assess this 
erroneous temperature excess in one-dimensional modeling of a mean spectra. 

In the present paper the usage of the H&K wings as model-construction 
diagnostics is calibrated in an effort similar to the analysis of Henriques and 
Kiselman (|2009|) and Henriques, Kiselman, and van Noort| (|2010j). The first 



goal of the paper is to estimate the reliability of one-dimensional modeling 
based on the Can H&K wings for characteristic fine structures of the solar 
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photosphere. The second goal is to quantify the temperature excesses produced 
by one-dimensional H & K wing modeling from unresolved spectra. 
Three different numerical simulations are employed here: 

• the 3D hydrodynamic (HD) simulation with the C0 5 BOLD code described 
by |Wedemeyer et aL|([2004]) . henceforth called 3DHD; 

• a 2D magnetohydrodynamic (MHD) simulation described bv lGadunl (|2000p 
and lGadun ~eTal\ POUT]) , henceforth 2DMHD; 

• a snapshot from a 3D magnetohydrodynamic simulation with the code of 
IStein and Nordlundl (|1998[) . henceforth 3DMHD. This snapshot was also 
used by ICarlsson efai\ ([2"UU2)l. Leenaarts et al. (j2"U06liI I2006b[) . Trit schler 
and Uitenbroek (|2006|), [ Uitenbroek, Tritschler, and Rimmele| (|200"T)) . and 
Henriques and Kiselman| p009[) . 

For each snapshot I computed emergent absolute intensities throughout H&K 
and then used these as quasi-observed H&K spectra to derive a best-fit one- 
dimensional photosphere model. I do not consider any actual H&K observation 
in this paper, but compare these simulation-result fits to the actual stratifications 
in the original simulation snapshots and compare the different simulations as 
well. 

In Section 2 I describe the simulation atmospheres, the spectral synthesis, and 
the procedure of one-dimensional model construction. The synthesized H&K 
spectra are presented in Section 3.1. The best- fit temperature stratifications of 
the average atmosphere inferred from spatially averaged synthetic H&K profiles 
and the evaluation of the temperature excesses are presented in Section 3.2. Best- 
fit temperature stratifications for granules, lanes and magnetic concentrations 
derived from spatially resolved H&K profiles are given in Section 3.3. Discussion 
and conclusions are given in Sections 4 and 5. 

2. Method 

2.1. Simulations 

The 3DHD simulation reproduces photosphcric granulation during a few hours of 
solar time. The simulation cubes measure (x, y, z) = 5600 x 5600 x 3110 km with 
an atmospheric extent of 1710 km in z above the zero "surface" level defined by 
mean optical depth logrs = at 5000 A. The grid step is 40 km in x and y; the 
height sampling in z starts at 46 km at the base and is 12 km for z > —270 km. 
I used 55 simulation cubes at 60 s time steps to obtain averaged H&K profiles. 

The 2DMHD simulation is representative of magnctogranulation with strong 
magnetic field, similar to network regions with fluxtubc field strengths reaching 
1.5-2 kG at the surface. The horizontal extent is about 4000 km, the vertical 
extent 735 km above and 1370 km below the zero "surface" level. The start-off 
initial magnetic field had bipolar geometry with loop connectivity and spatially 
averaged unsigned field strength 54 G, which increased with the simulation run 
time to a time-average of about 400 G at the surface. I used 52 snapshots at 
30-s intervals. 
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The 3DMHD simulation snapshot measures 6000 x 6000 x 5000 km with 
2500 km atmosphere height above logTs = 0. The grid step is 25 km horizontally, 
vertically 15 km in the upper layers and 35 km in the lower. The initial seed 
magnetic field was vertical and uniform with strength 250 G. The evolved snap- 
shot contains magnetogranulation with a few fluxtubes in intergranular lanes 
with field strengths of about 1700 G. Weak horizontal magnetic fields occur in 
the granules. 

2.2. Spectral synthesis 

I modified the SPANSAT code described by IGadun and Sheminoval ([1955)1 for 
LTE spectral synthesis of the overlapping H & K lines including superimposed 
blends. It computes intensity spectra containing any selection of overlapping 
spectral lines from a given one-dimensional model atmosphere. The spectral 
line parameters were taken from the VALD database ( |Kupka et al, 1999] ). The 
atomic parameters of the H & K Ca n lines are the same in Paper I, with calcium 
abundance A Ca = 6.38 and loggf = 0.134 for Can K, -0.18 for Can H. 

Micro- and macroturbulent velocities have small effect on the extended H & K 
wings. I did not apply micro and macroturbulence in the line synthesis from the 
simulations since their velocity field is given with the data cubes. In the sub- 
sequent fitting procedure I applied microturbulence in the line profile synthesis 
following IGurtovenko and Sheminoval (|1986p with the total velocity amplitude 
V t = 1.6 km s _1 determined as V t 2 = V£ ic + V^ ac . This approach simplifies the 
H & K profile calculation without loss of accuracy. 

Collisional damping by neutral hydrogen atoms was formerly a major uncer- 
tainty which has been resolved by the code of Barkle m and O'M ara (1998]) based 
on quantum-mechanical estimation. I used it here. 

In Paper I, which compared observed and computed H&K wing profiles, we 
had to cope with the "line haze" of unresolved weak lines in this part of the 
spectrum, for which we applied ad-hoc scaling factors 

/ = 1 + c(JV h /JVh-) (1) 

to the computed continuum opacity. The constants c were determined from com- 
paring synthetic far wings computed for disk center to the observed far wings. 
Since the present study does not involve actual observations this contribution is 
neglected here. 

In addition, no profile comparisons are made within 1 A from line center to 
avoid effects from coherent scattering. The effects of magnetic fields on H & K 
wing formation are neglected. All 2912 overlapping lines of neutral and single 
ionized atoms in the VALD list are entered as blends. In between these, I se- 
lected 28 relatively blend-free wing windows as the sampling wavelengths that 
define what I call "H & K profile" here. Such profiles were synthesized for all 
columns of all simulation snapshots. Then the average line profile was computed 
by averaging all these individual profiles per simulation. 

The major assumption in this computation is one-dimensional LTE radiative 
transfer along columns. IQwocki and Aue"rl(|1980p have shown that the unresolved 
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Mg k and Ca K wings are not sensitive to the horizontal transfer effects and may 
be used as potentially good atmospheric diagnostics in the multi-component 
(f .5D) sense without significant errors. More recently, ILeenaarts et al.\ (|20f Op 
have demonstrated that 3D radiative transfer effects only become noticeable at 
much larger height, as for example in the 3D evaluation of Na Di formation. 



2.3. One-dimensional model construction 



As in the classical H & K wing modeling I apply trial-and-crror fitting of the 
synthetic line profiles by changing the temperature stratification with height in 
the atmosphere as for a plane-parallel atmosphere. 

The initial start-off model gives the temperature T(h) and gas pressure P g (h) 
stratifications depending on the geometric height h (positive outwards). For 
the mean atmosphere modeling and for fitting granular fine structure I chose 
the HSRA-SP-M model as the start-off stratification. This is the HSRA model 
of Gingerich et al. (|1971[) modified by |Spruit| (| 1974[) and again modified in 



Paper I. Note that HSRA-SP-M contains very deep layers of the atmosphere. 
The PLA NEW reference model for a strong magnetic fluxtube by Solanki and 
Briglicvic (| 1992[) was chosen as the initial model for the modeling of small-scale 
magnetic elements. 

In each iteration the temperature stratification was changed manually at 14- 
17 geometrical height samples. Their distribution with height is irregular and 
was selected depending on the temperature gradients for the simulated object. 
The initial set of the height points (in kilometers) for the start-off HSRA-SP-M 
was the following: 864 (-6.7), 300, 200, 150, 100, 75, 50, 25, 0, -25, -50, -100, 
-150, -300, -409 (3.5). For PLANEW as start-off the set was 659 (-6.2), 300, 
200, 150, 100, 75, 50, 25, 0, -25, -50, -100, -150, -300, -306 (0.8). The highest 
and deepest optical depths are given between brackets, in logarithmic units. 
The temperature at other heights was computed by spline interpolation. The 
resulting best-fit stratifications extend down to the deepest geometrical depth. 
Since the deepest optical depth depends on the specific opacity K\(T, P g ) it 
depends on the simulated feature. 

At each temperature modification a new gas pressure stratification was com- 
puted assuming vertical hydrostatic equilibrium. The equation of hydrostatic 
equilibrium can be integrated in many ways. I used: 



P g (h) = Po exp [-1- / t^-dh'}, (2) 



RJh T{h< 

where g is the gravitational acceleration, R the gas constant, (J>o(h) the variable 
mean molecular weight taking into account the partial ionization of the various 
chemical elements. Po is defined as the gas pressure at the bottom of the atmo- 
sphere (ho) accounting for the adopted chemical composition. The SPANSAT 
code calculates the electron pressure P e as a function of depth from P g (h) and 
the equation of state for ideal gas including LTE evaluation of the ionization 
equilibrium for all pertinent elements according to the scheme by Mihalas (1967) 
with minor changes. In this work the chemical composition corresponds to data 
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Figure 1. (a) Mean H&K profiles synthesized from 3DHD (triangles), 2DMHD (asterisks), 
and 3DMHD (squares). The green observed profile with numerous blends is copied from the 
the disk-center Brault-Ncckcl atlas by Neckcl (1999). The intensity units were converted into 
brightness temperature, (b) Mean temperature stratification for the 3DHD (solid), 2DMHD 
(dashed), and 3DMHD (dotted) simulations. The spatial averaging was done over surfaces of 
equal optical depth rs . The temperature stratification of the quiet photosphere derived from un- 
resolved observed spectra HSRA-SP-M (Paper 1) and SRPM-305 (Fontenla, Balasubramaniam, 
and Harder, 2007} is presented by green and red curves, respectively. 



from the paper of Fontenla, Balasubramaniam, and H arder (2007 ). It is given as 
a table of abundance values for 104 chemical elements and it does not change dur- 
ing the fitting procedure. Thus, the free atmospheric parameters are (14—17) xT, 
Po in our one-dimensional modeling. 

As in other inversion methods each subsequent iteration minimizes the sum 
of the squared differences between the original and synthetic H&K profiles for 
28 wavelength points. The procedure of the minimization is complex due to the 
nonlinear dependence of the emergent intensity on the atmospheric parameters 
that control the radiative transfer. To succeed one must gain experience in how 
the intensity of the H&K profile responds to temperature changes of specific 
layers of the photosphere. 

Such fitting was done for the spatial average of each simulation and for single 
simulation pixels for which the modeling resolved specific fine structure (granule, 
intcrgranular lane, small-scale magnetic element). 



3. Results 

3.1. Synthesized spectra 

Figure QJi shows the averaged H&K profiles synthesized from all columns of 
the 3DHD, 2DMHD, and 3DMHD simulations. For the first two simulations 
the profile averaging is done over all columns of all snapshots, i.e., over both 
surface and time. In the third case the averaging is done over all columns of 
the single snapshot, i.e., only over surface. For comparison Figure [TJi also shows 
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Figure 2. Temperature against optical depth for all columns of a snapshot of the 3DHD and 
3DMHD simulations and for all columns of the time series of 2DMHD. The solid curve is the 
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Figure 3. Synthetic H&K profiles for all columns of all snapshots of 3DHD (a), 2DMHD 
(b), and a snapshot of 3DMHD (c). Temporal variations between profiles from all snapshots 
of the 3DHD (d) and 2DMHD (e) simulations. The symbols specify the mean value. In this 
and other figures the absolute intensity is given in watt cm -2 ster -1 A . 
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Figure 4. Best-fit H & K profiles (crosses) compared to the temporal and spatial averages 
of the synthetic profiles (rhombuses) from the 3DHD (a), 2DMHD (b), and 3DMHD (c) 
simulations. The lower panels of each subfigure show residuals AI = (/ — ^original) /^original; i- e -; 
the differences between the synthetic and original profiles in the relative units. The standard 
deviation a(I — / or iginal) are indicated in the lower-right corner of each subfigure. 

observed H & K profiles from the disk-center atlas taken by J.W. Brault with 
the NSO Fourier Transform Spectrometer and made available by INeckell (ff 999 ) 
after calibration to the absolute intensities of lNeckel and L abs (1984]). The three 
synthesized profiles agree reasonably with the atlas spectrum. However, the outer 
wings of the 2DMHD profile reach appreciably larger intensities than the other 
profiles. 

Figure |T|b shows the average temperature stratification for the three simula- 
tions, i.e., averaging both spatially and temporally for 3DHD and 2DMHD, 
spatially only 3DMHD. The spatial average is evaluated over surfaces with 
equal continuum optical depth T5 as in lUitenbroek and CriscuolTl (j2TjTTj) . They 
showed that such equal- T5 averaging represents emergent intensities better than 
averaging over planes of equal geometrical depth. 

As may be seen from Figure Q~b, in the deep layers below logrs=0 the mean 
gradient is smallest for 3DMHD and largest for 2DMHD. A steeper temperature 
gradient in deep layers produces larger emergent intensity in the outer wings. 
For comparison Figure [Tb adds the one-dimensional H & K line- fitting HSRA- 
SP-M model (Paper 1) used here to start the fitting procedure and the standard 
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Figure 5. Temperature stratifications of the best-fit models (dotted) compared to the 

equal- T5 averages of the simulations (red curve). The lower panels present the residual 
temperatures AT = T bcst _ flt - T original . 



one-dimensional c ontinuum-fitting SRPM-305 model of Fontenla, Balasubrama- 
niam, and Harder " i[2"HHT|r ? rheir temperature stratifications have lower gradients 
compared with the simulation averages. 

Figure [2] demonstrates the spread around the mean temperature stratifica- 
tion for each simulation. The spread caused by thermal inhomogeneities differs 
substantially between the three simulations. In high layers they are largest for 
2DMHD and remarkably small for 3DMHD; in deep layers they are more similar. 

Figures^, b,c show the spatial spread in the H & K intensities computed from 
the simulations, using all snapshots. In each case, the spread increases towards 
the wings with larger temperature spread at larger formation depth. 

Figures [3ji,e show the temporal variations between snapshots for 3DHD and 
2DMHD. The variability is largest for 2DMHD, probably reflecting growth in 
inhomogeneity because the magnetic field strength and oscillations increase with 
time. 



3.2. Best-fit stratifications from spatially averaged synthetic profiles 

This section presents best-fit results using the averaged H&K profiles. Figure H] 
shows the best-fit synthesized H&K profiles with the spatially averaged H&K 
profiles from all columns of 3DHD, 2DMHD, and 3DMHD, respectively. For each 
fitted profile the standard deviation is 



'original J 



\ 



l ™ 

-^(/(A,)- /original^)) 2 



(3) 



where the wavelength index i = 1, ...,n sample the wavelengths at which spec- 
trum is known, n = 28, /(A) is the intensity synthesized from the one-dimensional 
model, and /original (A) is the average intensity over the simulation. The largest 
values is a = 0.0008 for 2DMHD which has the largest temperature spread 
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Figure 6. Continuous intensity images at A = 3910 A for a 3DHD snapshot (a), the time series 
of 2DMHD (b), and a 3DMHD snapshot (c). The locations selected for feature modeling are 
marked with crosses and identifiers with designations: Gl and G2 for granule centers, GB1 and 
GB2 for granule borders, IG, IG1 and IG2 for intcrgranular lanes, LIG for a large intcrgranular 
lane, BP1 and BP2 for bright points, LMC1 and LMC2 for large magnetic concentrations, FT1 
and FT2 for centers of thin and moderate fiuxtubes. FT3 is a strong fiuxtubc. A cross in its 
center is called FT3-0 in the text; crosses to the left of center in the direction of the fiuxtubc 
periphery are called FT3-1, FT3-3, and FT3-5. The dotted line shows the selection of columns 
for the calculation of intensity and temperature fluctuations presented in Figure [T4l 



and gradient. In addition, the lower panels show residuals AI for quick as- 
sessment of the fit quality. These are normalized to /original per wavelength as 

AI (/ /original)//original- 

Figure [5] compares the corresponding best-fit temperature stratifications (dot- 
ted) to the simulation averages (solid red). The best agreement between the 
best-fit stratifications and the simulation averages occurs over \ogT$ = [0,-2]. 
Higher up the best-fit temperatures lie above the averages, especially for 2DMHD 
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Figure 7. Best-fit H&K profiles (crosses) compared to synthetic profiles from single 3DHD 
columns (rhombuses) for granule center (Gl, G2), granule borders (GB1, GB2), an intergran- 
ular lane (IG), and a large intcrgranular lane (LIG). The symbol and number coding arc the 
same as in Figure 3] 



due to the larger thermal inhomogeneity of 2DMHD (Figure Near logT5 « 
0.5 the best-fit stratifications lie above the averages for all three simulations. 
Figure [5] also shows that the temperature residual AT = Tbest-fit — ^original 
reaches a maximum of about 600 K for the 2DMHD stratifications, probably 
due to the growth of magnetic inhomogeneity. 

3.3. Best-fit temperature stratification from spatially resolved H&K profiles 

In this section a few columns corresponding to specific fine-structure features 
(granule center, granule border, intcrgranular lane, individual magnetic concen- 
tration) are selected for best-fit modeling. 

Figure |5] shows continuum images at A = 3910 A constructed from the three 
simulations with labels at the selected locations. The small-scale magnetic ele- 
ments in 2DMHD were selected with the help of vertical cross sections though the 
computational domain. In the 3DMHD simulation, they were selected according 
to the data of lCarlsson et al.\ (|2004p . 
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Figure 8. Best-fit temperature stratifications (dotted curves) derived from single 3DHD 
columns for the same locations as in Figure [7] Symbol and curve coding as in Figure [5] 



Figures [3 HI and [TT] show the best-fit results of the corresponding resolved 
emergent profiles. The best-fit synthesized H&K profiles (crosses) and the pro- 
files synthesized from the individual columns of the numerical simulation (rhom- 
buses) are in close agreement. The statistic estimations a (I — /original) m the 
lower-left corners and the residuals (in relative units) in the lower panels show 
the fit quality. 

Note the diversity between profiles from different types of granular structure. 
The steep wing slope becomes more gentle from a granule center to its borders 
and to a lane (Figure 0. The profile changes between the centre and periphery 
of strong fluxtubes are seen in Figures [Sg-j • The profiles from the bright points 
(Figures lllfc.f) are the narrowest. 

Figures [8j[T0j [12] show the best-fit temperature stratifications (dotted) and the 
corresponding original 3DHD, 2DMHD, 3DMHD columns (red). In general, the 
agreement is satisfactory. Some differences occur at temperature irregularities in 
the original stratifications (e.g., Figures ITUh.b.f.e) . As is seen from these figures, 
the temperature residuals AT are associated with the irregularities except for 
the large AT (200-400 K) at logr 5 = in Figure HDk.b. 
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Figure 9. Best-fit H&K profiles from 2DMHD for granule center (Gl, G2), intergranular 
lanes (IG1, IG2), thin fluxtubc (FT1), moderate fluxtube (FT2), and four locations in the 
strong fluxtube (FT3-0, FT3-1, FT3-3, and FT3-5). Symbol and curve coding as in Figured 



4. Discussion 



The discrepancies between the quiet-Sun atlas spectrum and the averaged profile 
from the 3DHD simulation without magnetic fields (Figure QJi) may have two 
causes. The first is the neglect of the "line haze" which affects the atlas profile but 
is not accounted for in the simulated profiles. Figure [T3] shows that inclusion of a 
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"line haze" in the H & K profile synthesis can significantly decrease the intensity 
in the H & K wings. Unfortunately, choosing an appropriate height-dependent 
fudge factor is not trivial. The second reason is that the temperature fluctuations 
in the simulation may not reproduce those in the real Sun. Evaluation of the 
simulation correctness is beyond this paper. 

The differences between the averaged profiles derived from the three simu- 
lations (Figure [TJl) can be due to the difference in the temperature inhomo- 
geneities, gradients, and magnetic field density (<|B|>). Figure [5b shows the 
largest temperature spread for 2DMHD which represents magnetogranulation 
with <|-B|> = 400 G in strong network. The large temperature spread in the 
higher layers of the 2DMHD increases the inner-wing intensities. It is caused by 
the presence of many magnetic network elements. Figure [2}) shows many sharp 
temperature drops and rises caused by the strongest magnetic fields within the 
fiuxtubes and at their hot walls at the peripheries. This follows from the analysis 
of vertical cross sections through the computational domain in the 2DMHD 
sequence. An example is Figure 10 of Paper 1. The temperature gradient of 
2DMHD is also the steepest one in Figure [Tb. 

The temperature fluctuations observed as granulation of the solar surface are 
created by penetrating convective flows. The thermodynamic properties of the 
granulation are described by IStein and Nordlundl (|1998p in detail. The temper- 
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Figure 11. Best-fit H&K profiles from 3DMHD for granule center (Gl, G2), intergranular 
lanes (IG1, IG2), bright points (BP1, BP2), and large magnetic concentrations (LMC1, LMC2). 
Symbol and curve coding as in Figure [4] 



ature and density of the plasma near the surface have approximately bimodal 
distributions. The intergranular lanes have low temperature, very low ionization, 
and high density. The granules have high temperature, high ionization, and low 
density. The granules occupy larger area than the lanes. Below the surface they 
occupy « | of the total area. Here the radiative cooling proceeds very rapidly. 
Since the H~ opacity is dominant and very temperature sensitive (w T 10 ), the 
steepest vertical temperature gradient results near the surface in the granules. 

The presence of magnetic field in a simulation also affects the granulation. 
IGadunl (|2000|) showed that the granulation pattern differs for magnetogranula- 
tion: the granules are larger and the shearing instability is stronger in the absence 
of magnetic fields. The field has a stabilizing effect - the granules are smaller 
and there are no pronounced horizontal shears. The distribution of granules 
in magnetogranulation differs from the nonmagnetic granulation by a greater 
contribution of bright features, which are associated with bright points in flux 
tubes. Therefore the granular contribution to the emergent intensity can also 
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Figure 12. Best-fit temperature stratifications from 3DMHD. Symbol and curve coding as in 
Figure [5] 
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Figure 13. Average H&K profiles computed without (triangles) and with (pluses) line-haze 
fudge factor of c = 0.14 in Eq. JTJ for the 3DHD simulation. The green observed profile from 
the Brault-Neckel atlas. 
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Figure 14. Fluctuations of the emergent intensity (<5/) in the far wing of the Ca II H line at 
A = 3910 A . The temperature fluctuations (5T) at the optical depth equal to unity inferred 
from a set of columns are signaled by the dotted line in a 3DHD snapshot (Figure [6ji) . 



differ between field-free granulation and magnetogranulation depending on the 
magnetic field density. 

In spatial averaging the temperature fluctuations have different effect on 
the temperature average and on the average emergent intensity, due to the 
nonlinea r dependence of the Planck function on temperature f Uitcnbrock and 
Criscuoli. l201ip . Indeed. Figure 1 141 shows intensity and temperature fluctuations 
(SI = (I - <!>)/<!> and ST = (T - <T>) / <T>) computed for a set columns 
through of a 3DHD snapshot (marked by a dotted line in the Figure H^i) ■ Here 
<T> is the mean temperature along the selected cut in this snapshot at op- 
tical depth equal to unity. <J> is the emergent intensity in the far wing of 
the Can H line at A = 3910 A averaged along the selected cut. The intensity 
fluctuations differ significantly in percentage from the temperature fluctuations. 
The relatively small temperature fluctuations ST = 5% produce large intensity 
fluctuations 51 = 50%. The intensity contributions to the average profile from the 
columns located within the granules are larger than from the columns located 
in intergranular lanes. This means that the averaged surface intensity from a 
simulation, i.e., averaged over the columns of a multi-dimensional simulation, is 
larger than the average intensity calculated from the averaged simulation. 

Figure [T5] shows the excess intensity between such averaged H & K profiles and 
the profiles synthesized from each spatially-averaged simulation. The continuum 
excess in relative units at A = 3910 A amounts to 7% for 3DHD, 38% for 
2DMHD, and 7% for 3DMHD. The intensity excess in the far wings is close 
to the continuum excess. Figure [3] demonstrates that the best-fit procedure 
leads to overcstimation of the temperature at continuum formation heights. 
The maximum temperature excesses are 194 K (3DHD), 619 K (2DMHD), and 
48 K (3DMHD). The larger the temperature inhomogeneities, the higher is the 
best-fit continuum and the larger the temperature excess of the best-fit average 
temperature stratification. Thus, in the deepest layers of the quiet photosphere 
the one-dimensional modeling of the mean atmosphere give overestimation of 
about 200 K, whereas in the strong network regions the overestimation can 
reach 600 K. 

The best-fit stratifications for the discrete fine-structure elements derived from 
the synthetic H&K spectra in Figures [81 ITUI [T2l recover the original temperature 



SOLA: pap.tex; 10 July 2012; 0:20; p. 17 



V.A. Shcminova 



0.6 



I" 0.4 

c 

£ 

0.2 




0.6 

If 0.4 
c 

0.2 




0.6 



i 0.4 
c 
o 
c 

0.2 


3920 3940 3960 3980 
wavelength [A] 

Figure 15. Averaged line profiles synthesized from the 3DHD, 2DMHD, 3DMHD simulations 
(dotted curves) and the profiles synthesized from their average temperature stratifications 
(solid), as well as the corresponding continua. 

stratifications satisfactorily. Notice that the fitting procedure applied here recov- 
ers successfully both a temperature stratification which decreases monotonically 
with height and stratifications with a mild or steep temperature rise in higher 
layers (Figures l5F.HUe.h.f.i.l51i.fTUb.fT^b.e.f.g.hV The question arises whether it is 
possible to identify a temperature rise from the H & K profiles alone. Figures [7J 
l^lllllshow that H & K profiles calculated for temperature-rise stratifications have 
higher intensity near the line core the for a stratification with decreasing tem- 
perature. Probably, this diagnostic can be employed indicator of an atmospheric 
temperature rise. 

In Figures 151 HUlfT^l thc temperature residuals AT = Tbcst-fit — ^original reflect 
the intrinsic uncertainties of the method used. For example, Figure [TOb shows 
large (about 400 K) deviations from the original models at logrs = 0, i.e., 
large deficiencies of the method. Tests showed that adding radiation pressure, 
magnetic pressure, and turbulent pressure in Eq. ([2j) makes no significant dif- 
ference. Since the chemical composition is the same for the profile synthesis 
and the best-fit modeling it does not affect the results. The real cause for such 
deviations may be a effect of thermal irregularities in the vertical direction on the 
best-fit stratification. The effect is tested in Figure [TBI The smooth temperature 
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Figure 16. Temperature (a) and gas pressure (b) stratifications for a granule (Figure HOb l 
as well as contribution functions (c) to the emergent intensity at A = 3910 A and line profiles 
(d): original (red curves), smoothed (black), and best fit (dotted). Contribution functions are 
given in logarithmic units. 



stratification (black curve) differs from the original (red curve) at locations where 
there are sudden temperature changes with depth. Such thermal irregularities 
are not recovered by the method used. In Figure [TBb the gas pressure P g (h) 
(black) derived from the smooth T(h) with Eq. ^ has excesses compared to 
the original one (red) at the locations corresponding to sharp temperature gra- 
dients. The synthesized H&K profile (black) with the smooth T(h) therefore 
docs not match the original profile (red) in Figure ITSH. The corresponding black 
contribution function for the intensity at A = 3910 A (Figure [16b ) is slightly 
lower than the red one due to increase of the damping (7 w d w ~ Pg) and the line 
opacity k(T,P 3 ). It causes lower emergent intensity in the far wings. To obtain 
a best line-profile fit the temperature increases and the pressure decreases (see 
dotted curves). Therefore the best-fit temperature (dotted curve) demonstrates 
the excess about 400 K in Figure [TBk. 

Complex non-monotonic temperature stratification such as in Figure fTUH .e.f 
also cannot be reproduced by the manual fitting procedure used here. In addition, 
in the manual fitting the initial T(h) and especially P g (h) stratifications in the 
deepest layers must be as close as possible to the real ones. So, for example, 
in the case of magnetic elements it is not possible to obtain best-fit profiles 
with the HSRA-SP-M model as initial guess, whereas the PLANEW fluxtube 
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model gives good results. These are shortcomings of the manual technique which 
diminish wh en using automatic inversion codes (e.g., Bellot Rubio, Ruiz Cobo, 
and Collados, 1997; Frutig er et q<.[l2000p that search for a best fit by means of a 
Marquardt nonlinear least-squares algorithm. They can locate the minimum of 
the merit function more precisely than a manual procedure. The retrieved tem- 
perature stratification by the automatic inversion codes is largely independent of 
the choice of the initial guess atmosphere. Note, it is a limitation for the manual 
technique, but not for the H&K wings photospheric diagnostics. Therefore, it 
is desirable to develop such inversion codes for fitting H&K profiles. 



5. Conclusion 

The extended wings of Ca II H & K supply relatively easy diagnostics to extract 
quantities that describe vertical stratification of the solar photosphere. I have 
calibrated such one-dimensional modeling using three numerical simulations of 
solar fine structure, containing granular convection with different amounts of 
magnetic flux. I obtained best-fit temperature stratifications from synthetic 
H&K wing profiles computed from these simulations. Quantitative evaluation 
permits the following conclusions. 

The different effect of thermal inhomogeneities on the average temperature 
stratification and on the average emergent intensity due to the non-linear Planck 
function sensitivity leads to overestimation of temperatures derived in such 
one-dimensional modeling. For quiet-Sun conditions this overestimation reaches 
about 200 K in the continuum formation layers. The larger the inhomogeneity, 
the larger the excess. The largest overestimation occurs for one-dimensional 
modeling at field strengths similar to the 2DMHD simulation, and can be as 
large as 600 K for strong magnetic fields when average profiles are analyzed. 
Such temperature overestimation may also occur in the upper photosphere due 
to large temperature contrasts between hot fluxtubc walls and cool fluxtube 
insides. 

The H&K wing diagnostics provide better recovery of the actual temperature 
stratifications within discrete fine structure elements such as a granule center, 
an intergranular lane, a bright point, a fluxtube, a large magnetic concentration. 
Here errors occur only where there are sharp vertical gradients and non-smooth 
irregularities in the original temperature stratification. 

The range of effective formation of the emergent intensity in the wings of Ca II 
H&K extends from logrs = —2 to 0.2. Beyond this depth range the reliability 
of the H&K wing diagnostics is less. 

In conclusion, one-dimensional H&K wing modeling using high-resolution 
data is a reliable technique to construct best-fit models for small-scale features 
in the solar photosphere. 
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